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Observations of the superconducting energy gap in 
the far-infrared spectra of conventional superconduc- 
tors have played an important role in establishing the 
Bardeen-Cooper-Schrieffer (BCS) theory of superconduc- 
tivity 0, and provided an evidence for importance of 
strong-coupling effects predicted by Eliashberg and Hol- 
stein theories 

Despite numerous studies |^ Q IE 111 > for Nb the influ- 
ence of the strong coupling on optical properties is not as 
clear as it is for Pb or Hg, for instance. Most recent re- 
ports on experimental results of the surface impedance in 
the microwave region are rather contradictory. Anlage et 
al. 0] and Klein et al. described their data basically 
by the BCS theory, while Marsiglio et al. were not able 
to fit their results either by the simple BCS model or by 
calculations according to the Eliashberg theory citeeliash 
with the behavior of a'^F{uj) evaluated from tunnelling 
measurements 

Only a few experimental studies 0, 0| determine 
both components of the complex conductivity, cr{uj) = 
ai{uj) + i(j2{w), as a function of frequency in the most in- 
teresting frequency region - around the superconducting 
energy gap 2 A, that allows to make quantitative com- 
parison with theoretical strong-coupling calculations. 

In one of such experiments 0, both, the real cti(ijj) 
and the imaginary a2(L^) parts of the complex conductiv- 
ity have been directly determined from the transmission 
and the phase shift spectra of high-quality thin niobium 
films on transparent substrates. The measurements have 
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been performed at several temperatures above and be- 
low Tc in the frequency range 5 to 30 cm~^, precisely the 
region where the superconducting gap lies. It has been 
shown that the overall frequency dependence of the con- 
ductivity can be described using the BCS formalism and 
assuming finite scattering effects. At the lowest temper- 
atures, however, some deviations from the BCS predic- 
tions have been observed in the frequency behavior of the 
complex conductivity. The deviations in the real part of 
a{uj) could not be explained by the strong-coupling ef- 
fects. In a subsequent comment [lH . an attempt has 
been undertaken to explain the reported deviations by 
the granular structure of the sputtered Nb films. These 
films have been considered as a two-component system 
with an intrinsic BCS type conductivity ctbcs and a 
grain-boundary conductivity abn of Josephson type. 

In this communication we show that the experimen- 
tally measured in Ref. [loj data for niobium, can be very 
accurately described using the first-principles approach 
developed in Ref. No Josephson junctions inside 

the sample are needed to be considered. 

This first-principles approach (for a review see 
allows to calculate the electron and the phonon spectra of 
metals and the spectral functions of the electron-phonon 
interaction. As a result, many related to the electron- 
phonon interaction properties, such as the temperature- 
dependent electrical and thermal resistivity, the plasma 
frequencies, the electron-phonon coupling constants, etc., 
have been computed in Ref. [l^ . 

For Nb, it has been shown that the electron-phonon 
contribution to the dc resistivity p{T) is well described 
by the first-principles calculated values of the plasma 
frequency ujpi—d.Al eV and of the transport electron- 
phonon coupling constant Atr=1.17 [l^ . 



The resistivity p{T) of niobium can be written as: 

p{T) = i^7(r), (1) 



where 7(T) is the relaxation rate 



The first term here is the relaxation rate due to the tem- 
perature independent impurity scattering. The second 
term describes the temperature dependent contribution 
from the electron-phonon interaction. 

At considerably high temperatures (T > 0.29d, here 
9_D is the Debay temperature) Eq. El can be simphfied 
to: 



7(r) = 7i,„p + 27rAfrT, 



(3) 



where Xtr is the transport constant of the electron- 
phonon coupling: 



A. 



(4) 



Fig. n shows temperature dependent part of the dc 
resistivity of the thin magnetron-sputtered Nb film llR l 
in comparison with the literature data for bulk Nb 
The residual impurity resistivity of the film is po = 
4.1 X 10~^r2cm [lOj. As it is seen from the figure, the 
experimental data obtained on a very thin (150 A) film 
do not demonstrate significate difference from the bulk 
results. The slope of the film curve is somewhat stepper 
than this for the bulk curve. A little decrease of a plasma 
frequency (by no more than 10% of the bulk value) can 
describe this difference. 

The critical temperature of the film, Tc — 8.31 K ,10i|, is 
also somewhat below the bulk value (9.3 K), although for 
films of this thickness, it is rather high The influence 
of the film thickness on the properties of thin films in the 
normal and superconducting states is certainly beyond 
the scope of this work. For detail discussions on this 
topic we refer to Ref. |0| and references therein. Here we 
just emphasize that the deviations of all the parameters 
for the film measured in Ref. ^3 from the bulk values 
are within ±10%, thus the investigated film can pretty 
accurately represent the bulk properties. 

From the theoretical point of view, these small differ- 
ences between the film and bulk properties can be under- 
stood, for instance, in terms of a simple approach devel- 
oped by Testardi and Matteiss This approach in- 
troduces smearing of the first-principles calculated band 
structure of metals due to the lifetime effects, which are 
certainly very essential in thin films, and allows to calcu- 
late the density of states A^(0) and the plasma frequency 
as function of the phonon (pph) and the residual (po) 
resistivity. Although the Testardi-Matteiss approach is 
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FIG. 1: Temperature dependence of phonon part of dc resis- 
tivity for thin-film [l^ and bulk [l^ niobium. 



oversimplified, it gives the right direction of the change 
in the film properties compared to the bulk. 

For calculation of the frequency dependent conductiv- 
ity of Nb, we have used the formulas, obtained firstly by 
Nam 20] . Here we present them in a shape convenient 
for numerical computations: 



cr(w) 



} 



^ 1 + n{uj — uj')n{uj') — a{uj — Lu')a{uj') 

e{uj - LU') + e{uj') + 2ijirnp 

1 — n{LL! — uj')n*{Lj') + a{uj — Lj')a*{Lj') 
£{uj - cj') - e*{uj') + 2i7„„p 

,1 — n{u! + w')n(w') — a{u! + uj')a{uj') 
e{uj + uj') + e{uj') + 2i^imp 
1 + n[uj + w')n*(w') + a(w + u}')a*{u}') 
£{uj + Lo') - e*{uj') 4- 2i7„„p 



} 



x_Re 



dio' tanh( '^ ^ ) 
1 — n{LO + Lo')n{u}') — a{u} -\- Lu')a{uj') 

e{uj + Lu') + e{uj') + 2i"/irnp 

' , UJ - Lu' 

du tanh( ) 
„ 1 — niuj — Lu')n*(Lu') + aiuj — uj')a*{Lo') . 

^ - - L 

£{U0 ~ Uj') - £*{uj') + 2i^irnp 

where / is the Fermi distribution and 



(5) 
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FIG. 2: Frequency dependence of the real (upper panel) and 
imaginary (bottom panel) parts of the complex conductivity 
of Nb, as calculated from the first principles. Compare to Fig. 
4 of Ref. flo(|. As an example, we replot the experimental 
(72 (<^) data from this figure at T = 6 K. 



FIG. 3: Examples of the transmission and phase shift spectra 
of a 150-A-thick Nb film on sapphire substrate. Open symbols 

- spectra, calculated from the first principles. Solid symbols 

- experimental data from Ref. [lo| . 



a{uj) 



^cj2 _ A2(w) 
A 



(6) 



The renormalization function Z{lu) and the energy gap 
A(a;) have been calculated from the Ehashberg equation 
0, which has a well-known form, and we do not repro- 
duce it here. In our calculations, wc used the Eliash- 
berg spectral function a'^F{uj) obtained from the first- 
principles calculations in Ref. [l^ . 

In order to match to the experimentally measured crit- 
ical temperature of the film, we had to slightly increase 
the Coulomb pseudopotential jj,* in the Eliashberg equa- 
tion, in comparison with the bulk value. As films are 
usually more disordered than bulk samples, and as the 
disorder is known to increase the value of fi* , this increas- 
ing seems to be well grounded ^1]. Due to the reasons 
discussed above, while computing the Eq. |31 the plasma 
frequency has also been slightly changed from its bulk 
value of 9.47 eV to 9.2 eV. 

The real and imaginary parts of the complex conduc- 
tivity, as obtained from our calculations, are shown in 
Fig. 12 for several temperatures below Tc- At each tem- 
perature, the calculations have been made for 20 - 30 
frequency points in the interval from 5 to 30 cm~^. As 



temperature goes down, the energy gap clearly develops 
in the (7i(a;) spectra, (72 (ti^) demonstrating a characteris- 
tic l/tj-divergency at — > 0. 

To compare our computations with the raw experimen- 
tal data, we calculate the transmission and the phase 
shift spectra of a Nb film on sapphire substrate (as has 
been measured experimentally in Ref. ^^I)- F'or these 
calculations, the optical parameters of the substrate and 
the thicknesses of the film and of the substrate have been 
taken from the same reference. The calculations were 
performed at all temperatures where the experimental 
data were available. An example of these calculations 
together with the experimental data are shown in Fig. |2| 
for T = 6 K. Clearly, the computed spectra beautifully 
describe row experimental results. 

In Fig. [3 as large open symbols, we show (T2 data, cal- 
culated in Ref. ^IQ] from the experimentally measured 
transmission and phase shift spectra for T = 6 K. Our 
first-principles calculations fit these experimental points 
very well, while the weak-coupling calculations have been 
shown to significantly deviate from the experimental val- 
ues, especially at the lowest frequencies (Fig. 5 of Ref. 
^3)- Thus, the effect of strong coupling is clearly seen 
in the frequency-dependent conductivity in Nb. 

We should note, that at the lowest temperatures and 
low frequencies, the real part of the complex conductiv- 
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ity, calculated in Ref. ^3 from the measured transmis- 
sion and phase shift spectra, does not coincide with our 
first-principles calculations. We suggest this deviation 
is caused by the difficulties in solving the inverse task 
- calculation of ai{uj) and a2{oj) from transmission and 
phase shift. At the lowest temperatures and at frequen- 
cies below the gap, the imaginary part of conductivity is 
much large than the real part. Consequently, the exper- 
imentally measured dynamic response (transmission and 
phase shift in our case) is mainly determined by 172 ('^), 
while (Ti (w) has quite a minor influence on experimental 
spectra. That causes large error bars in determination of 
ai{u)) from experimental data. 

Since the relaxation rate (which can be determined 
from Eqs. ^-13 is much larger than the gap value, or, 
in other words, the mean free path is smaller than the 
coherence length, the film is in the dirty limit and the 
following expression for the penetration depth is valid: 



A"^(T) = cr„ • 27r/TO 



A2(w)tanh(w/2r)dw 
w2 - A2(a;) ' 



where tT„ is the normal state conductivity: 



47r7i 



(7) 



(8) 



imp 



For T = 4K, that gives A — 880 A, in a good agreement 
with the experimental result, A^^^* = 900 A^. 

To make better comparison with the experimental re- 
sults of Ref. 01 , we now compute the London penetra- 
tion depth. By the London penetration depth we mean 
the penetration depth that would be observed in our su- 
perconductor, if it were in the clean limit. For this goal, 
we shall put jimp = in Eqs. [SJ what gives: 



27r 



Re 



A2(w)tanh(w/2T)dw 
(tj2 _ A^{Lo) f/^ ■ ReZ{Lu) 



(9) 



For T = 4 K we obtain Al = 354 A, again in a per- 
fect agreement with the experimental result, A^^'' = 350 

A[i3. 

In conclusion, we have demonstrated that the row ex- 
perimental data on the frequency dependent conductivity 
of a thin Nb film, obtained in Ref. can be accu- 

rately described by the first-principles calculations made 
for bulk Nb with just a small (±10%) varying of such pa- 
rameters as the plasma frequency and the Coulomb pseu- 
dopotential. These deviations are most likely connected 
to the disordered nature of films, and can be understood 
in terms of the a simple theoretical approach No 
Josephson junctions inside the sample are needed to be 
considered to describe the experimental findings. For the 
first time, in our knowledge, the first-principles calcula- 
tions of the complex conductivity of a superconductor, 
have been successfully compared with the experimentally 



measured data in the frequency region around the super- 
conducting gap. 
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